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Abstract 

The observed flavor structure of the standard model arises naturahy in "spht fermion" models 
which localize fermions at different places in an extra dimension. It has, until now, been assumed 
that the bulk masses for such fermions can be chosen to be flavor diagonal simultaneously at every 
point in the extra dimension, with all the flavor violation coming from the Yukawa couplings to the 
Higgs. We consider the more natural possibility in which the bulk masses cannot be simultaneously 
diagonalized, that is, that they are twisted in flavor space. We show that, in general, this does not 
disturb the natural generation of hierarchies in the flavor parameters. Moreover, it is conceivable 
that all the flavor mixing and CP-violation in the standard model may come only from twisting, 
with the flve-dimensional Yukawa couplings taken to be universal. 
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I. INTRODUCTION 



One of the motivations to extend the Standard Model (SM) is to explain the fermion 
flavor structures. It is likely that there is a more fundamental theory that produces the 
observed masses and mixing angles in a natural way, namely, without small dimensionless 
parameters. One such framework uses split fermions to generate the small numbers j3| • The 
basic idea is to localize the SM fermion fields at different locations in compact extra dimen- 
sions. Then, the four dimensional (4D) Yukawa couplings between left handed and right 
handed 4D fermion fields are exponentially suppressed by the overlap of the corresponding 
zero mode wavefunctions. In general, such a split fermions setup induces small and hier- 
archical 4D Yukawa couplings without imposing any extra symmetries. In addition, one 
can account for proton stability by separating quarks and leptons in the extra dimension. 
Specific realizations, phenomenological implications and experimental s igna tures of the split 
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fermions framework can be found in 

Fermion localization works as follows [l|, 12, ISl- Consider, for simplicity, a model with 
one infinite extra dimension. (For the more realistic case of a finite extra dimension see, for 
example, .) The model contains one bulk scalar, the localizer, which is assumed to 

get a vacuum expectation value (vev) which depends on the extra dimension coordinate, z. 
Thus, the five dimensional (5D) Dirac spinors have two mass terms: a ^-independent bare 
mass and a z-dependent mass term due to the couplings to the localizer. For \Efj, a generic 
Dirac field (with i = 1, 2, 3 as the generation index) and $(-2), the localizer vev, these mass 
terms read 

£ = ¥Mi,{z)^\ Mij{z) = rrii, + K.^z) , (1) 

where rriij and Xij are z-independent Hermitian matrices. 

In the past is was always assumed that rriij and Xij can be diagonalized simultaneously. In 
that case the problem of obtaining the 4D observables is significantly simpler. We can carry 
out the Kaluza-Klein (KK) decomposition of the 5D fields in the basis in which the mass 
matrix is diagonal. Each zero mode, which is interpreted as the corresponding SM chiral 
fermion field, then has a wavefunction in the extra dimension fi{z). These wavefunctions 
must satisfy the condition derived from the 5D Dirac equation: 



- Mi{z)\ fi{z) = 0, Mi = [mu + \ii^{z)] 
where = d/dz. The solution is 



fi{z) = /i(0)exp 



r MAz')dz 
Jo 



(2) 



(3) 



where /j(0) is an arbitrary vector that specifies the values of the solutions at z = 0. For 
example, in the case where the localizer vev is linear, the zero mode wavefunctions are 
Gaussians P, 0|. These Gaussians peak at the points where Mi(z) = 0. The fact that 
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the wavefunctions which correspond to different generations are locahzed at different points 
leads to the required small overlaps. 

In this work we relax the assumption that the parameters in the Lagrangian should 
be aligned in flavor space, that is, we assume that rriij and Xij cannot be diagonalized 
simultaneously. We call the unaligned case "twisted", and the aligned case "untwisted". We 
will focus on answering the following two questions related to the attractive features of the 
split fermions framework: 

(i) In the presence of twisting, can one naturally suppress operators which involve fields 
in different representations? This is required to account for the proton longevity. 

(ii) In the presence of twisting, are the hierarchies between operators which involve fields 
from the same representation still natural? This is required to account for the flavor 
puzzle. 

As we demonstrate below the answer to both of the above questions is positive: The presence 
of twisting does not spoil the basic appealing features of the split fermions framework. 

Once we understand how twisting affects hierarchies, we can ask whether twisting may 
also be useful for model building. Again the answer is positive. One example is related to 
the fact that twisting induces new CP violating sources. This was used in js| to construct 
a new type of leptogenesis model. Below we demonstrate how the standard model flavor 
mixing and CP violation may arise purely from twisting. Another possible application is 
related to a solution of the strong CP problem, which will be discussed in a separate work. 

II. THE TWIST - BASIC FORMALISM 

We consider the most general Lagrangian of Eq. (Q). We would like to find the profile of 
the zero modes and calculate the 4D observables in terms of the 5D Lagrangian parameters. 
The answer is nontrivial because the effective mass matrix of the fermions Mij{z) is a z- 
dependent Hermitian matrix. There is no global SU{3) flavor transformation which brings 
Mij{z) to a diagonal form simultaneously at every z. Mij{z) can be diagonalized by a z- 
dependent special unitary matrix U{z), however this does not leave the 5D kinetic terms 
invariant. This is to be compared to the untwisted case where U is ^-independent. Formally 
we can introduce a local measure of the twist by 



When Rij{z) is non-zero, a twist is present in that region. Only when Rij{z) = for all 
values of z does the general (twisted) case reduce to the flavor- aligned (untwisted) case. 
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(4) 
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In the general twisted case the KK decomposition for a vector (in flavor space) fermion 
field is given by 



Z) = Y.UI)U^)PMX,) + llUR)^c.{z)Pnrc.{x,) 



(5) 




where are the known four dimensions, {fi^ia and {fji)ia are ^;-dependent wavef unctions, 
and Pl {Pr) is the left (right) handed chirality projection operator. Our notation is such 
that 5D (4D) fields are denoted by capital (lowercase) letters. We can think about the Latin 
(Greek) indices as labeling the flavor space in 5D (4D). Note that the 5D wavef unctions, /]^ 
and /£, are generic z-dependent matrices in flavor space, as compared to the untwisted case 
0], where they were diagonal and hence simply functions. 

We are interested in the zero modes since they correspond to the SM fermions. Consider, 
for example, the wavefunctions of the left handed zero modes, fja{z). These must satisfy 
the condition derived from the 5D Dirac equation: 



where j = 1 . . . 3 stand for the three components of a single wave function, while a = 1 ... 3 
labels the three different solutions to this equation and thus corresponds to the SM flavor 
indices. The solution may be written formally in a straightforward way: 



where P stands for the path ordered product and /ja(0) is an arbitrary 3x3 matrix that 
specifies the values of the different solutions at 2; = 0. In principle, we can choose /ja(0) 
such that the three vectors fn, fi2 and fi^ constitute a set of orthonormal eigenvectors: 



In practice, we can just choose the vectors to be linearly independent and then get an 
orthonormal set by applying the Gram-Schmidt procedure. 

The nontrivial difference between the twisted ((Tj) and untwisted © case is that with 
twisting one cannot factorize the solution of the zero modes into a form of flavor-space times 
z-space. At each point, each solution is a vector in flavor space, but the twisting forces it 
to rotate (twist) as it moves along the extra-dimension. This follows directly from the fact 
that with twisting the matrix Mij cannot be diagonalized simultaneously at all z. 

Unfortunately, in the most general twisted case, equation ((Tj) is not very instructive 
(although it can be used for numerical computations). In certain cases even with twisting 
we can find explicit solutions. For example, in two generations, the two coupled first-order 
differential equations © can be combined into a single second order equation. Then, if M 
has a simple enough form, we may be able to find solutions. For example, if M depends 
linearly on z, the solutions are Kummer functions (see Appendix^. In three generations 
the composite equation is third order and is in general unsolvable, even with a linear M. 



5ijd^ - Mij{z) fja{z) = 0, 



(6) 




(7) 




(8) 
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In order to obtain the 4D observables we start from the 5D couphngs of the fermions to 
the SM Higgs field. Consider for example the couplings of the quarks doublets, Qi, to the 
down type quark singlets, Dj, 

S = J (fxYf^HQ.Dj , (9) 

where H is the SM Higgs field and the Dirac structure is suppressed. The 5D Yukawa 
couplings, Fj^, are assumed to be arbitrary parameters without any specific flavor structure. 
Performing the KK reduction, assuming that the profile of the Higgs vev is flat, and keeping 
only the zero modes, we obtain the standard 4D action 

S = j d^xyiphq^dp. (10) 

The dimensionless 4D Yukawa couplings are given by 

yU = J dzY.j f!:{z) ff,{z) , (11) 

where (f^) is the wavefunction of the left (right) handed quark doublet (down type 
singlet) and the sum over i and j is implicit. 

For simplicity, here and in it what follows, we work with rescaled parameters. That is, 
we scale constants and wavefunctions to the appropriate power of the fundamental scale to 
make them dimensionless. 



III. GENERAL PROPERTIES OF THE ZERO MODES 

In the following two sections we demonstrate that twisting does not destroy the essential 
appealing features of the split fermions framework, namely, localization and separation of 
the fermion zero modes. 

First, consider localization. We confine our discussion to the case were the eigenvalues 
of Mij(z) are monotonically decreasing functions of z (as is the case in the models of 3])- 
In the untwisted case the peaks of the zero modes are located at the points where one of 
the flavors has a zero bulk mass. This implies that we can separate fermions in different 
representations, for example quarks and leptons, to forbid proton decay. In this section we 
show that a similar localization holds in the twisted well. 

To see this, we show that since the eigenvalues of Mij{z) are monotonic functions of z, 
we can define a "localization region" . We define this region as the region between the two 
points Za and Zb such that for z < Za {z > Zb) all the eigenvalues are positive (negative). 
The magnitude of each of the zero modes, l/ap = X)i flafia: decays outside the localization 
region as is shown below. 

It is sufficient to study the norm of the wavefunction because the overlap between the 
norms provides an upper bound on the overlap between the actual wavefunctions. 



dzn:{z)f-p{z) < / dz\r^{z)\\f;^{z)i (u) 
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which serves as a bound on the effective 4D couphngs. Using eq. ^ we obtain 

d,\Uz)\' = 2Y^m,{z)\fU^)\'. (13) 

k 

where there is no sum on the index a and mk{z) are the eigenvalues of Mij{z). We denote by 
"^max(^) ["^min(2;)] the maximal [minimal] eigenvalue of M at z. Then, we can place bounds 
on the right hand side of equation (fT^ 



<E^^(^)I/-WI' <^max|/.(^)P. (14) 



We see that the norms of the zero modes fall at large positive and negative values. In 
particular, 



l/.(^)r<exp 



dzlrrir^J^z!) 



\faiZb)\'' i0TZ>Zh. (15) 



Note that mmax < over all of the integration range. A similar bound may be placed in the 
z < z„ region 



l/.(^)l'<exp 



rza 

-2 / dzmrn\n{z) 

J z 



\UZa)\^ hiZKZa, (16) 



where m^^a > all over the integration range. We see that the norm of the wave function 
indeed decays outside the locahzation region. 

The above arguments hold separately for each SM representation. Thus, different rep- 
resentations can be separated if their localization regions do not overlap. In models where 
the localizer vev is linear, like that of [l], mmax and mmin are roughly linear far away from 
the localization region. In that case the bounds in equations fll5j) and (fTHIl imply that far 
from the localization region the norms of the zero modes are suppressed exponentially (as 
Gaussians). The generalization to the more realistic models in which a stable scalar configu- 
ration is not a monotonic function of z is straightforward. In these cases one can divide the 
extra dimension into regions where the scalar is monotonic and apply the above analysis for 
each of the regions separately. Thus we have answered question (i) from the introduction: 
fermions in different representations can be naturally separated. 



IV. THE ADIABATIC APPROXIMATION 

Though we have shown that exponentially small overlaps are easily achieved between dif- 
ferent representations, in order to solve the flavor puzzle hierarchies within a representation 
are required. To see whether this is the case in our framework one needs a better handle on 
the profiles of fermions zero modes. As we already mentioned, the general solution to the 
zero mode wavefunction equation is not very useful in this respect. Here we show that in 
many cases one can make an approximation in which the governing physics is clear. 

The zero mode equation ^ bears resemblance to the time dependent Schrodinger equa- 
tion, idtlip) = Hlip), with it replaced by z. The z dependent M corresponds to a time varying 
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Hamiltonian. One of the useful approximation methods for solving a time varying Hamilto- 
nian in quantum mechanics is the adiabatic approximation where the wavefunction is 
assumed to be an instantaneous energy eigenstate at all times. 

It is useful to make a similar approximation in our case. The assumption we make is that 
the independent solutions of eq. (0) are each governed by a single eigenvalue of the mass 
matrix at each point in the extra dimension. We thus expect each solution to be aligned 
with one of the eigenvectors of M{z) at each point. 

We begin by writing an ansatz for the solution to the equation which satisfies the guideline 
mentioned above. The ansatz for the left handed zero mode profile is 



fta (z) = ^ exp 



dz' rria {z') 



V,^ (z) , (17) 



where Via{z) are the z dependent normalized eigenvectors of the twisted bulk mass matrix 

Mij{z)Vjc, (z) = m^{z)Via{z) with V*,Vip = 6^^ . (18) 

Within this approximation each zero mode is localized around the zero of a single {z de- 
pendent) eigenvalue of the mass matrix. The adiabatic ansatz is thus a straightforward 
generalization of the untwisted solution. In both cases the wavefunction of each profile is 
governed by a single function which is an eigenvalue of the bulk mass matrix and therefore 
in both cases the wavefunctions are localized. In the adiabatic limit the only added feature 
introduced with twisting is that the wavefunction points in different directions in flavor space 
for various values of z. 

To what extent is the adiabatic approximation valid in generic models? The standard 
condition in quantum mechanics literature is that the approximation holds so long as the 
quantity 

{V^\d,M\V2) ^ (19) 
(mi — ' 

is small throughout the evolution of the system jly]. Our case however, is somewhat more 
subtle due to the fact that the evolution of the wavefunctions with 'time' is not unitary. 
The rate at which the true solution deviates from the adiabatic one at a certain point is 
indeed proportional to the quantity in eq. (fT^ but is also depends on the values of the 
wavefunctions at that point. 

We did not try here to fully formulate and derive the necessary conditions for the validity 
of the above approximation. One can, however, get good intuition by numerically comparing 
it with the exact solution, eq. (|7j). This allows us to learn about the accuracy of the adiabatic 
approximation. As an example consider the following mass matrix for the two generation 
case 

MW^-^r' (20) 



where we work with rescaled parameters such that z and M are dimensionless. In the upper 
part of figure ^ we plot the two eigenvalues of the mass matrix and in the lower part of the 
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FIG. 1: The two solid blue lines in the upper figure are the eigenvalues of the mass matrix M{z) 
given in eq. (|2fl)) . The localization region defined in the previous section lies between the two 
vertical dotted lines. The exact solution to the norm of the wavefunctions is in the lower figure 
in solid lines while the adiabatic approximation is in thick dashed green curves. Note that the 
departure of the adiabatic solution from the exact one occurs where the two eigenvalues approach 
one another. 

figure [U we plot the norm of adiabatic solution that is induced by the eigenvalues (green 
dashed curves). For comparison the adiabatic curves are plotted on top of the exact solution 
(black solid lines). 

We see that in the case presented in figure [T] the adiabatic approximation is very good. 
It noticeably deviates from the exact solution only in the region where the two eigenvalues 
approach one another as one would expect from eq. (fT^ . This departure only occurs once 
both profiles have decayed (or have started to do so). Consequently, physical quantities are 
hardly affected and the approximation holds to a very good accuracy. 

We expect this behavior to be generic in models that have well localized and separated 
wavefunctions. This is because separation implies that at the points Zi, where any of the mass 
eigenvalues vanishes, the mass differences are large. Thus, non-adiabatic contributions are 
supressed near Zj. If the eigenvalues are changing smoothly the non-adiabaticity generically 
occurs well outside the localization region. We have indeed observed that increasing the 
separation (and hierarchy) between flavors improves the approximation. We have checked 
numerically several cases that lead to hierarchical Yukawa matrices and found that the 
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adiabatic approximation works very well in these cases. 

V. HIERARCHIES AND MIXING FROM TWISTING 

So far we have shown that even with twisting, localization and separation within a rep- 
resentation naturally occur. In this section we discuss twisting as the only source of flavor 
mixing and CP violation. We consider below the following two questions: 

(a) Can all the flavor mixing and CP violation come from twisting? That is, can we take 
the 5D Yukawa matrices to be universal, so that all of the SM flavor structure comes 
from twisting? 

(b) Can the above scenario be realized naturally? 

We show that the answer to the above two questions is positive. We first demonstrate 
numerically that twisting can serve as the only source of flavor mixing. Then we support this 
observation, for more generic cases, using the adiabatic approximation. Finally, we present 
a toy model in which this situation occurs naturally. 

Let us begin with the numerical example. We denote the mass matrices as 

Ml, {z) = - Til^z . (21) 

with r = Q,U, D. We assume that the 5D Yukawa couplings are proportional to the unit 
matrix. Then, the Lagrangian of the model is given schematically by 

C = Y^SijHQi Uj + Y^'SijHQi D, + M^QiQj + M^^^UiUj + MgDiDj . (22) 

When and rj^j can be diagonalized simultaneously this model is untwisted. In the most 
generic case, however, when no symmetry is imposed, the mass matrix is expected to be 
twisted. For our explicit example we choose the following parameters for the bulk masses 

M«(^) 
M^(^) 

where we work with rescaled parameters such that z and M'^ are dimensionless. We obtain 
hierarchical fermion masses 

^ = 2.1 X 10-\ ^ = 8.0 X 10-^ (24) 



3Az-A 2.3 \ 
2.3 2.6^-6/ ' 



4.1^ + 1 1 \ 

1 3.1^-3/ ' 

'4.5^ + 4 1 \ 

1 2.9^ + 2 J ' 



(23) 
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and finite mixing 

sinec7 = 0.19. (25) 

Note that this is only an example. We did not try to search for a mass matrix which precisely 
generates the observed 4D flavor parameters. We use this example only to demonstrate the 
fact that 4D flavor mixing and hierarchical masses can originate only from the twisting. 

While we gave an explicit example only for one simple case, the conclusion holds much 
more generally. To see this, we use the adiabatic approximation. Using Eqs. (jllj) and (jl7p 
we see that 



z 



V?^ {z) {z) . (26) 



The exponential factors in ()26|) teach us that we expect a hierarchical structure for the 
4D Yukawa couplings. In addition, for the realistic case of three generations, the product 
V^Vj^* is just a product of two SU (3) matrices. This shows that flavor mixing is induced. 
Furthermore, these SU{3) matrices generically contain 0(1) phases. Thus, ^ckm, the CP 
violating phase in the 4D CKM matrix, is expected to be of order unity as observed. 

The above numerical study exemplifies the possibility that all of the SM flavor conversion 
is achieved due to the twist in the bulk, and not due to the 5D Yukawa couplings. Below we 
construct a toy model which naturally realizes this idea. Consider a model for quarks with a 
non-Abelian horizontal flavor symmetry SU{3) p on a 5D orbifold M4 x S1/Z2. The fermions, 
Q^,W, are fundamentals of the flavor group where all other SM fields are singlets. The 
Higgs field, the left handed component of and the right handed component of and 

are even under the Z2 orbifold symmetry, and all the other fields are odd. Because of 
the orbifold symmetry, a bulk mass for the fermions is forbidden. But an effective mass can 
be generated from the vev of Z2 odd bulk scalars that are SM singlets. For our example, 
we include an adjoint (octet) $ij and a singlet (p of SU{?>)f- The Yukawa couplings to the 
Higgs, on the other hand, are allowed and are proportional to unit matrices due to the flavor 
symmetry. 

Generically, there is a potential for the scalars in the bulk. This potential naively gen- 
erates an untwisted scalar vev since varying vevs are usually not the lowest energy config- 
uration. Thus, we also include boundary terms that change this naive expectation. When 
the symmetry is explicitly broken on the boundaries, the competition between the bulk and 
boundary terms can force the vev of $ to be twisted. 

To be explicit, consider the case in which the boundaries preserve only an SU{2) subgroup 
of the bulk SU{3)p. Under this SU{2), decomposes into a scalar <l> = $33, a fundamental 
$i = ($13, $23); an antifundamental and an adjoint with l,m = 1,2. For each of 
the above fields we can write a boundary term 

_^brane^ ^ {d^^f + bl {d^^lf + cl{d,^ ij'] + . . . (27) 

B=0,ttR 
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where higher order (stabihzing) terms are omitted for simphcity. The parameters ao,7r_R, 6o,7r_R 
and Co,7r_R are real and generically different on the two different branes. The terms above 
cause non-zero vevs to be developed for the derivative of $, $i and at 2; = 



where Him is an order one Hermitian matrix and we used an SU (2) transformation to bring 
the vev of dz^i to its special form. Similarly, at z = nR the derivatives get non- vanishing 
vevs, which are generically different from those at 2; = 0. Consequently, the mass matrix at 
both boundaries are not aligned and a twist is generated in the bulk. 

We assume that v\l <^ 1 where vp is the typical vev of the bulk scalars and A is the 
effective cutoff of the 5D theory. Then, to leading order in Vp/A"^, the 5D Lagrangian is 
given by eq. with 



where c\ and C2 are unknown constants. We see that the 5D Yukawa couplings are propor- 
tional to the unit matrix and that M[j cannot be globally diagonalized at each z because it 
is twisted. Due to the fact that the bulk scalars are odd under the Z2 orbifold symmetry 
the correction to the universal Yukawa matrices is suppressed by 0{vp/A'^). There are addi- 
tional higher dimension brane Yukawa terms which are volume suppressed and are therefore 
negligible. We conclude that in this toy model the 4D flavor violation is, to leading order, 
only due to twisting. 

VI. DISCUSSIONS AND CONCLUSIONS 

We found that in many ways neglecting the twist is a simplifying assumption. Twisting 
does not change the fact that split fermions generically produce hierarchical Yukawa cou- 
plings. At first neglecting the twist seems unjustified: it changes the symmetry breaking 
pattern of the theory. For example, suppose there is a single fermion representation in the 
bulk, with 3 flavors. When M{z) oc I3 there is an enhanced U{3) flavor symmetry in the 
theory. Including a diagonal M{z) breaks the U{3) symmetry down to U{1)^. Including a 
generic twisted M{z) breaks it further to U{1) In untwisted models this last step occurs 
only due to the standard model Yukawa interactions. 

In the untwisted case the fact that we get hierarchical Yukawa couplings and small mixing 
can be understood from symmetry considerations as follows. The U{1)^ — > U{1) breaking 
is due to non-diagonal 5D Yukawa couplings, which connect fields of different SM represen- 
tations, and therefore, the last stage of symmetry breaking occurs between objects that are 
separated in the extra dimension. Thus, the 4D Yukawa couplings are small since they are 
suppressed by the small overlap of the separated zero modes. That is, the U{lY remains 
as an approximate symmetry, and it is restored once the zero mode wavefunctions are far 
away from each other. 





(28) 



(29) 
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In the twisted case the fact that an approximate f/(l)^ symmetry is obtained in the low 
energy effective theory is less obvious. Nevertheless, we claim that this remain the case due 
to the same reason, namely, symmetry breaking occurs between objects that are separated 
in the extra dimension. In other words, even with the twist, we have shown that the zero 
modes are localized and separated. This implies that the 4D Yukawa couplings are small 
since they are proportional to the small overlap of the different zero modes. Just like in the 
untwisted case, when the separation is very large, the U{lY symmetry is restored. Thus, 
the effect of the twist is only to add new sources of flavor mixing and CP violation. 

In conclusion, the main result from our study is that the twist does not affect the gen- 
eral appealing features of the split fermions framework, that is, the possibility of naturally 
creating hierarchies without symmetries. Furthermore, it opens a possibility in which the 
observed CKM mixing and CP violation arise from twisting and not from the 5D Yukawa 
couplings. 
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APPENDIX A: LINEAR POTENTIAL: EXPLICIT SOLUTION 

There is one case where we could find an explicit analytic solution to the zero mode 
wavef unctions. This is the case of a two generation model with an infinite extra dimension 
and a localizer with a linear vev. Here we only sketch the derivation [17], and discuss some 
of the properties of the solution. 

In the case under study the mass matrix defined in (0) can be written as 

Mij{z) = rriij - r]ijZ . (Al) 
It is convenient to chose a basis in which niij is real, rjij is diagonal and rjn > 7722. Then, 
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eq. © is given explicitly by 



/i(^) + (^112; - mu)fi{z) - 77112/2(2) = 0, 
/2(^) + (^222; - m22)f2iz) - mufiiz) = 0, 



(A2) 
(A3) 



where the prime denotes a derivative with respect to z and the index a was dropped. The 
solution is 



/i(z) = p{z) [c,mut{z)M{a, 3/2, w{z)) + c^M^a - 1/2, 1/2, w{z))] 
h{z) = p{z) [c2m,2t{z)M{a + 1/2, 3/2, w{z)) + ciM(a, 1/2, w{z))] 



(A4) 



where M is the Kummer function and 



p{z) 



exp 



(mil - rjii 



27] 
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Vu -r]22 + m 



12 



t(z) 
ti'(2;) 

The normalization condition is 



2(?7ii - 7722) 

"^11 - "^22 

z , 

Vu - V22 

^-^^^tizf. 



dz = 1. 



(A5) 



(A6) 



The set of equations ()A2|) and ()A3|) has two independent solutions related to the two inde- 
pendent constants Ci and C2 that appear in ()A4|) . This degree of freedom in the solution 
corresponds to the index a of fia- One can choose the integration constants such that the 
two wavefunctions are orthogonal 

[dzj2M^)Mz) = 0. (AT) 

In practice we ensure the orthogonality by applying the Gram-Schmidt procedure on a pair 
of non-orthogonal wavefunctions. 

Arriving to the solution in ()A4|) is straightforward. We first use flA2|) and plug it into 
()A3|) to arrive at a second order differential equation for /i. Using fi{z) = g{z)p{z)t{z) this 
equation is 



which is the Kummer equation [l8|. The general solution to ()A8|) is 



g{w) = Qa M{a, 3/2, w) + gb 



+ {3/2-w)-f^-ag = 0, 

w 

solution 

M(a- 1/2, 1/2, w) 



(A8) 



(A9) 
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where ga and gb are two independent constants. To get the solution for /2 we used the 
following properties of M jist - 



w 



7 



dM(a, 7, w) a , ^/ , , x 

\' " > = -M a + 1,7 + 1, w , 

aw 7 

M(a + 1, 7 + 1, w) = M(a + 1, 7, w) - M(a, 7, w). 



aM(a + 1, 7 + 1, w) = (a — 7)M(q;, 7 + 1, w) + 7M(a, 7, if). 



(AlO) 

(All) 
(A12) 



One can check that in the mi2 limit the solution of the twisted case ()A4|) reduces 
to the solution of the untwisted case. In that limit a ^ 1/2 and we recall the following 
properties of the Kummer function 



M(0,c, h) 
for arbitrary c and h. Then we get 

(mil - ^711 zf 



1, 



M(c, c, 6) 



/i(^) 



C2 exp 



27] 



11 



/2W 



Ci exp 



(m22 - r]22 z) 



2V: 



22 



(A13) 



(AM) 



Taking the two independent solution to be those where either ci or C2 vanish, we get the 
two Gaussian solutions of the untwisted case 
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